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Abstract 

A broader perspective is suggested for the study of higher dimensional 



cosmological models. 



1 Introduction 



■ — i ■ 

' In this context one might imagine "the view from above" to refer to the vantage 



point of the extra dimensions of spacetime somewhere "up there" above us, as 
in the fiber bundle diagrams one often sees in discussions of higher dimensional 
theories. However, in using this phrase here I have in mind something quite 
different. 

No matter how many extra dimensions we consider for spacetime, they all 
live in a 2-dimensional world of paper, blackboards, computer screens, overhead 
transparencies, ... In fact an incredible quantity of 2-dimensional space has 
been filled with discussions of these higher dimensional models. It is often easier 
to keep grinding out the next small step in this production line, but sometimes 
one has to step back and take a more global look at the situation, to rise up 
above this 2-dimensional detail (in our 3-dimensional world) and think about 
broader questions. How do things fit together? Are there properties which are 
general for the many special cases considered? What properties of these special 
cases will continue to have validity in a more general setting? 



1 



In these brief remarks I cannot say many things which deserve being said. 
Being somewhat of an outsider in the field, I must also be careful not to say 
things which reveal my own ignorance of and lack of familiarity with much of 
the current work in higher dimensional models. Having recently studied some 
properties of classical symmetric cosmological models in higher dimensions, 1 
this invited talk pushed me to try to make some sense of some of the details 
scattered in the literature. 

Since time has been short, I cannot explain exactly how everything fits to- 
gether, but I hope to convey a general picture which if pursued may put many 
of these results into perspective. In so doing I would suggest that perhaps a 
little more attention should be paid to how specific calculations fit into the 
larger scheme of things. Certainly it is much more satisfying when many little 
things fit together as particular cases of a single bigger thing. This idea charac- 
terizes present day theoretical physics and the importance it gives to unifying 
symmetries. 

Unfortunately the search for exact solutions of 4-dimensional gravitational 
theories has more often than not focused on special cases without putting them 
into perspective. In higher dimensions there is much more room to play games. I 
think the situation demands some restraint. Some effort should be made toward 
"looking at the forest rather than the trees." The game of looking for exact or 
even qualitative solutions of gravitational field equations all too often finds itself 
at the tree level when an aerial view of the forest is really what is needed. 

Unification of the fundamental forces has led people to adopt the Kaluza- 
Klcin idea 2-5 of formulating theories on higher dimensional spacetimes. At first 
the extra dimensions were introduced as a mathematical device for obtaining 
dimensionally reduced unified theories on 4-dimensional spacetime. These addi- 
tional dimensions, tied to symmetries of the spacetime, allowed the collapse of 
the theory to an equivalent 4-dimensional one. Eventually ideas progressed to 
the point where the extra dimensions were taken seriously. The "fiber" symme- 
try of these fiber bundle formulations of 4-dimensional theories was abandoned 
and such theories became truly higher dimensional ones. ' 

This immediately created a problem. If the fundamental theory of mat- 
ter is really higher dimensional in such a nontrivial way, one must explain the 
present effective 4-dimensional spacetime and the apparently 4-dimcnsional the- 
ory which describes nature. Some mechanism is required to collapse the addi- 
tional dimensions to a point where they are not observable. This could be ac- 
complished if they are associated with spatially compact dimensions with very 
small circumferences, but one needs some way for this "compactification" to be 
a natural consequence of the theory itself. 

An interesting possibility for this compactification mechanism is the dynam- 
ical behavior of the gravitational field itself within these theories. The present 
4-dimensional phase of the universe is very anisotropic, given the necessarily 
large difference in length scales associated with the "ordinary" and "extra" di- 
mensions. Chodos and Detweiler 8 first suggested this anisotropy might be a 
consequence of the evolution of the gravitational field. Examining purely grav- 
itational 5-dimensional theories with flat homogeneous space sections charac- 
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terized entirely by "free" anisotropic expansion, 9 they noted that classical cos- 
mological solutions are allowed in which preferential expansion of three spatial 
directions occurs. However, such solutions are not selected by the theory. One 
needs matter fields to obtain such a selection effect. 10 Such fields may be added 
in ordinary Kaluza-Klein theories, while supersymmetry dictates what matter 
fields exist in locally supersymmetric theories (supergravity, superstrings). 11-19 

With this in mind classical solutions of the Bose sector of higher dimensional 
theories have been studied to throw some light on the nature of the ground 
state of such theories. Such studies fall into two categories, characterized by 
the behavior of the gravitational field. In the first, solutions are considered in 
which the spacetime is the product of two static symmetric component manifolds 
representing the ordinary spacetime and the internal spaces, the discussion being 
more particle physics in flavor (e.g., reference 13). In the second, evolving 
models are considered, typically from the viewpoint of a relativist somewhat 
less concerned with the particle physics motivation for considering such models. 

As a classical relativist I wish to address the second kind of investigation in 
the spirit of Freund, 10 Dcmaret et al 20 ' 21 and others. 22-27 Many particular exact 
or qualitative cosmological solutions of various higher dimensional theories have 
been discussed. Although only for very special conditions can one obtain closed 
form solutions, such solutions may describe a single isolated dynamical effect 
of the full equations which when accompanied by other degrees of freedom no 
longer allows exact solutions but perhaps exhibits a similar qualitative behavior, 
or at least one which can be understood in the context of the coupled behavior. 
Properties exhibited by special solutions are probably only interesting to the 
extent that they relate to the general case in this way. 

Certain exact solutions in fact exhibit a "typical behavior" during certain 
phases of the evolution. Such solutions can be understood in terms of additional 
properties related to a qualitative formulation of the field equations. From a 
purely mathematical viewpoint, the existence of closed form exact solutions 
often is accompanied by underlying properties of the differential equations that 
are rarely understood in first obtaining them and which help to relate various 
particular solutions to each other. Understanding how exact solutions fit into 
a bigger picture is perhaps more important than any particular exact solution 
found by chance or "cleverness". 

These ideas are also useful in understanding the idea of "gravitational chaos" 
both in the general cosmological context and in the context of the homoge- 
neous models. 28-31,23-26 This question is studied as a purely gravitational phe- 
nomenon relevant to the classical initial singularity independent of particle the- 
ory questions. Although quantum effects at the Planck length probably make 
the question irrelevant, it is an interesting mathematical problem. 

Limitations of time and space permit only a hint of a more global point of 
view illustrated by a few key examples. Hopefully such a discussion, necessarily 
abbreviated, can still be valuable in establishing the importance of such a point 
of view for this branch of cosmology. 



4 



2 Spatially Homogeneous Models 



The kinds of cosmological solutions considered in higher dimensional theories in 
connection with possible properties of the "ground state" all share one feature: 
they are spatially homogeneous. 32 The (i-dimensional spacetime contains a 1- 
parameter family of (d — l)-dimensional spacelike hypersurfaces (let D = d — 1) 
which are orbits of an isometry group G. The space sections are therefore 
"homogeneous spaces" diffeomorphic to the coset space G/H , where H is the 
isotropy group at some point of spacetime, namely the subgroup of isometries 
which leave that point fixed. 

When H is the identity subgroup or a discrete invariant (i.e., normal) sub- 
group, "space" is isometric to a group manifold (namely the group G/H) with 
a Riemannian metric invariant under either the left or right translations of the 
group into itself, the choice being a matter of convention (left in this article) . If 
H is discrete but not invariant, this is still true locally. All of these cases may 
be referred to as the simply transitive case, where G acts as a D-dimensional 
translation group of a Z?-dimensional Riemannian manifold. 

In the multiply transitive case, H is continuous and G acts as a group of 
translations and certain rotations of the space sections, each of which is iso- 
metric to the left coset space G/H with a Riemannian metric invariant under 
the natural left action of G. If H is an invariant subgroup, G/H is a group 
and hence this occurs as a special case of the simply transitive case for which 
additional symmetry exists; the case where H is not invariant will be referred 
to as nontrivial. 

In either case the metric on a given spatial section is completely determined 
by its value at a single point and hence the field equations reduce to ordinary 
differential equations in the parameter labeling the spatial slices, easily inter- 
preted as a time function on the spacetime. (This assumes that the matter fields 
either share the symmetry of the metric or at least have enough symmetry to 
be compatible with the field equations.) By choosing the coordinate lines or- 
thogonal to the spatial sections, the spacetime metric can be expressed at least 
locally in the form 33 

^g = -~N(t) 2 dt®dt + g ab (t)uj a ® uj b (2.1) 

on the manifold R x (G/H), where {uj a } are 1-forms on G/H dual to the frame 
{e a } a =i,...,D characterized by structure functions C a bc = uj a ([e a ,eb\) on G/H. 
These are constants in the simply transitive case where {e a } is a left invariant 
frame (i.e., a basis of the Lie algebra g of G), but functions on G/H in the 
nontrivial multiply transitive case. In the first case the matrix g = (g a b) is an 
arbitrary positive definite symmetric matrix, while in the second case it must 
satisfy additional linear constraints. 18 ' 19 Similar constraints in the first case may 
lead to additional spatial symmetries. 

The choice of lapse function N(t) determines the parametrization of the 
family of homogeneous spatial sections. The proper time r, defined by dr = 
N(t)dt = uj 1 - corresponds to unit lapse. In expanding models with an initial 
singularity, r is usually chosen so that this singularity occurs at r = 0. Defining 
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BIANCHI TYPES 


nonscmisimple 


semisimple 




noncompact 


compact 




abelian 


nonabelian 


class A (unimodular) 


I 


II, VI , VIIo 


VIII 


IX 


class B (nonumimodular) 




III, IV, V, VI fc , VII,, 





Table 1: The isomorphism classes of 3-dimensional Lie groups/ algebras. The 
compact/noncompact division refers to the simply connected covering groups of 
each type, while unimodularity refers to the matrix representation of the adjoint 
group of each type. 



lu° = dt leads to the spacetimc 1-forms u> a dual to the frame {e a } a= o,i,. ..,£>, 
where eo = d/dt. 

In practice one is usually interested in spatial geometries which decompose 
into an orthogonal product of an ordinary 3-dimensional homogeneous space 
and a compact internal homogeneous space 

G/H - M 3 x (G mt /H rnt ) . (2.2) 

M3 is a Bianchi manifold 34 in the simply transitive (group manifold) case and 
the Kantowski-Sachs manifold 35 in the nontrivial multiply transitive case, where 
M 3 is of the form G/H with G ~ Gi x SO(3, R), Gi the abelian group Rot Si, 
and H ~ SO(2, R). (One might refer to M 3 as a BKS space.) A classification of 
symmetry types in the Bianchi case amounts to a classification of 3-dimensional 
Lie groups/ algebras; 34,36 such a classification is summarized in Table 1. 



3 Field Equations 

Since we wish to focus on the gravitational properties of higher dimensional 
models, we consider the Einstein equations with all other fields lumped into 
the energy-momentum tensor. This assumes that the gravitational part of the 
classical field equations can be represented in this way, thus excluding "higher 
derivative" theories. These equations 

= M afj ee iWffWGetf - T aP ) , \ {d) g\ 1/2 = Ng 1 / 2 (3.1) 

(where metric quantities without the leading superscript refer to the spatial 
metric) naturally split into evolution equations = M a b which evolve the spatial 
metric and constraints on the solutions of those equations 

= H ee 2N~ 1 M J -± = 2M°i (super-Hamiltonian constraint) 
= H a = 2N~ 1 M ± a ee 2M° a (supermomentum constraint) . 
The Einstein equations may also be written in "Ricci form" 
= P a p ee ^R a0 - E aP , E aP ee T aP - (d - 2)- 1 T\ g aP , (3.3) 
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a form very convenient for static solutions in discussions of spontaneous com- 
pactification. However, for dynamical solutions a compromise is preferable, 
namely the Ricci evolution equations = P ab , supplemented by the Einstein 
constraints rather than more equations involving second time derivatives. 

To understand the usefulness of this compromise for spatially homogeneous 
solutions, it is convenient to introduce the Hamiltonian formulation. The Ein- 
stein evolution equations written in first order form are the Hamiltonian equa- 
tions for the spatial metric derived from the Hamiltonian H = NH with the 
lapse treated as an independent variable; here Tr ab is the momentum canonically 
conjugate to the spatial metric g ab 

= 7r ab + d/dg ab (NH) - NQ ab 

= jr ab + NdH/dg ab - NQ ab + d/dg ab (\nN) = M ab + Hd/dg ab (\nN) . 

(3.4) 

(The term Q ab arises in the case of nonunimodular symmetry groups where the 
spatial Einstein tensor force is no longer derivable from the scalar curvature 
potential but has a nonpotential component or from similar problems with the 
source energy-momentum. 37 ' 38 ) If N is treated as independent of g ab , the fi- 
nal term vanishes and one obtains the Einstein evolution equations. Letting N 
depend explicitly on g ab leads to evolution equations which differ by a term pro- 
portional to the super-Hamiltonian constraint. The choice of unit lapse N = 1 
made in proper time time gauge leads to the Einstein evolution equations, while 
the choice N = {const)g x l 2 made in Taub's time gauge 32 (Misner supertime 39 
time gauge) leads to a set of equations which are equivalent to the Ricci evolu- 
tion equations. 

Why is the Taub time gauge useful, i.e., why is it useful to choose the lapse 
proportional to the spatial volume factor? For the answer, one must examine 
the differential geometry of the Einstein equations in Hamiltonian form. The 
metric configuration space is the space M of positive definite symmetric matrices 
g = {gab) or a submanifold of this space. The Hamiltonian 

H = NH = NT + NU (3.5) 

consists of a kinetic term associated with the metric variables and a potential 
energy term composed of a spatial curvature term and a source term 

U = U curv + U source , -dU + Q ab dg ab = -g 1/2 (G ab - T ab )dg ab , 

U curv = -g 1/2 R , -dU curv + Q curv ab dg ab = -g 1 ' 2 G ab dg ab , 

fjsource = _ 2 . g 1 / 2 T ± ± , _ dI jsource + Qsource ab^ = g^l^ T ab dg ab . 

(3.6) 

The exterior derivative d is the one on M. , while both the source and the curva- 
ture forces may involve nonpotential components. (One may also simply treat 
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the source force entirely as an external force.) The kinetic energy term 

NT = \N- x Q ahcd g ah g cd (velocity phase space) , 

= NQ~^ cd -K ab ir cd (momentum phase space) , (3.7) 
TT ab = d/dg a b(NT) (Legendre transformation) , 

is just the one associated with the DeWitt 40 metric Q on M, apart from the 
rescaling by the lapse function 

N~ 1 Q = JV-y^aOy^ - g ab g cd )dg ab ® dg cd . (3.8) 

Allowing the lapse to have a factor depending explicitly on g a b is equivalent 
to conformally rescaling the DeWitt metric by that factor. 39 In particular the 
Taub time gauge choice eliminates the spatial volume factor from the DeWitt 
metric, thus simplifying the kinetic terms in the evolution equations as much 
as possible and also eliminating a term involving the total kinetic energy which 
appears in the Einstein form of those equations due to the derivative of the 
spatial volume factor. This explains the advantage of using the Ricci evolution 
equations accompanied by the Taub time gauge. 

Neglecting the potential energy equivalent to restricting one's attention to 
the vacuum abclian case where the spatial curvature vanishes, the evolution 
equations describe geodesies of the DeWitt metric Q in proper time time gauge 
and geodesies of the rescaled DeWitt metric g~ x ^Q in Taub time gauge. The 
Hamiltonian constraint in this case requires that the kinetic energy vanish and 
therefore be a null vector. (In fact the null geodesies are conformally invariant.) 

This immediately raises the question of the signature of the DeWitt metric. 
It is a Lorentz metric which endows M. with a natural splitting into space 
and time which itself plays an important role in the geometry of the Einstein 
equations. Let 

g = e 2 «g, g = dctg = e 2Da 1 detg=l (3.9) 

decompose g into the "scale invariant" conformal metric g and the conformal 
scale factor e a . Dynamically in the spatially homogeneous models e a describes 
the overall expansion or contraction of the spatial sections while g describes the 
anisotropy. The DeWitt metric in these new variables is explicitly Lorentzian 
with a playing the role of a natural time variable 

g- 1/2 G = -£)(D -l)da(g) da + Tr(g- 1 dg(g)g- 1 dg) . (3.10) 

The second term is Ricmannian since at g — 1 it leads to the trace of the 
square of a symmetric matrix and such an inner product is positive definite on 
the space of symmetric matrices. 

The diagonal submanifold A4d C M. is in fact a flat Lorentz spacetime with 
respect to the rescaled DeWitt metric. At diagonal points g£> one may define 
= \ lngu and introduce a basis {e^} a=\,...,d-i of tracefree diagonal matrices 
normalized by the condition Treves = D(D — 1)6 ab, and let /3 = /3 A eA, in 
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which case {a, f3 A } are orthonormal cartesian coordinates with respect to the 
suitably rescaled DcWitt metric 

[D(D - l).g 1/2 ] _1 aU D = -da ®da + S AB dp A ® d/3 B . (3.11) 

The vacuum abelian diagonal dynamics in Taub time gauge has as solutions 
straight null lines in a-0 space (unit a-speed trajectories in space) with these 
orthonormal coordinates linear in the Taub time t. (It is convenient to choose 
N = 2D(D — l)^ 1 / 2 to obtain the usual factor of one half in the kinetic en- 
ergy.) The diagonal components of g are therefore exponential in the Taub 
time. Transforming to proper time time gauge leads to power law dependence 
on the proper time r 

g D = e A *+ B = g D(0) diag(r 2 ^, . . . ,r 2 ^) , (3.12) 

the latter form recognizable as the famous Kasner solution, the Kasncr expo- 
nents satisfying the well known identities 

The Ricci evolution equations in Taub time gauge for diagonal metrics are 
explicitly 

(lng)- = -[2D(D - l)] 2 • 2g(R - E) (3.14) 

where R = (R a b) and E = (E a b). When the righthand side is absent, one has 
free motion with the metric exponential in the Taub time. This remains true 
for those degrees of freedom along which the force term on the righthand side 
has no component. 

As an aside, it is worth mentioning how spatial curvature and nonzero source 
energy density change the causal character of the solution curves. The super- 
Hamiltonian constraint T = —U determines the square of the velocity (or mo- 
mentum) of the system. In the case of positive source energy density, U may be 
negative and the motion spacelike only for sufficiently positive spatial curvature. 
In a contracting phase this can lead to a "bounce" rather than a cosmological 
singularity, while in an expanding phase it can lead to a point of maximum 
expansion and recollapse. Such "turnaround phases" are not allowed when the 
spatial curvature is negative and the source energy density is positive, where U 
is positive and the motion always timclike, with g or equivalently a monotoni- 
cally increasing (expansion) or decreasing (contraction), respectively starting or 
ending at the "frontier" g = (a — ► -co), where the spatial metric is singular. 40 
The intersection of a trajectory with the frontier may or may not represent a 
physical cosmological singularity. The variable a G (-co, oo) is a natural time 
variable for M ; during a phase of monotonic expansion or contraction, one can 
choose it as the time variable t on the underlying spacetime, a gauge introduced 
by Misner 41 (a-time time gauge). 
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4 Exact Solutions 



A useful example to illustrate some of these ideas and others are the Taub 
solutions 32 in d = 4 spacetime dimensions. This class of metrics is defined by 

C\, = C 2 3i = n« , C\ 2 = , gD = e 2a e 2 3+e ++ /J-e_) ) 

e+=diag(l,l,-2), e_ = V3diag(l,-1,0) , 0~ = 5(n<U)fa , 

(4.1) 

and their dynamics in Taub time gauge is described by the Hamiltonian (ex- 
pressed in velocity phase space) 

H = ±(-a 2 + $+ 2 + /T 2 ) + 6n^ 2 e 4 ^- 2 ^ - 24n^ n« e" 2 ^ ^ . (4.2) 

The Lorentz transformation 39 

(a, /?+) = 3- 1/2 (2a -(3+,-a + 2(3+) (4.3) 

leads to a completely decoupled Hamiltonian 

H = i(-S 2 +t +2 + t 2 ) + 6n( 3 > 2 e- 4 ^ + - 2An^n^e 2 ^ 
= -Ha + Hp + + Hp- = 



(4.4) 



consisting of three 1-dimensional scattering problems with exponential poten- 
tials and constant energies, restricted only by the constraint that the appropri- 
ately signed sum of the individual energies vanish. The solution of an exponen- 
tial scattering problem leads to hyperbolic or trigonometric functions 

H x = \{x 2 + = E X - t = J dx{2E x - iie^Y 1 ' 2 , 

or defining 7 = (2E x y/ 2 , (4.5) 

\vx _ f V 1 ' 2 !' 1 cosh(± 7 i4) n > , E x > 
ln^-y- 1 wah(^\u\t) M<0, E x € 5i . 



e 2 



(a,(3 + ,P ) are inertial coordinates of the rest frame of the n^ 3 -* 2 potential, 
which moves with speed d(3 + /da = \ in /3 space. For the Bianchi type IX 
case, characterized by n^rS 1 ^ > 0, the hyperbolic cosine solution is relevant, 
interpolating between the asymptotic free positive and negative exponential 
solutions (for the metric components) at t = ±00; the unit velocities in /3 space 
of the asymptotic solutions are related by a simple reflection in the rest frame 
of this potential. Letting — > contracts the group to Bianchi type II, 
eliminating the "tachyonic" n^n^ potential (it moves with speed d(3 + /da = 2 
in space) and allows free motion parallel to the 2 potential (the a and 
(3~ directions). For the Bianchi type VIII case, characterized by n^n^ < 0, 
only the positive energy solutions are relevant for the a motion (the potential 
is negative) due to the super-Hamiltonian constraint. 
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evolution 
equation (set =0) 


spatial 
curvature 


2nd time 
derivative 


quadratic 
decoupling 


quadratic 
elimination 


linear 
choice 


A R A + *G ± ± 


n W 


A 


n = -1 


m = 1 


(1,-1) 


Rb 


n (3)2 


B 


to = 2 


n = 1 


(2,1) 


4 R b - 2( 4 G ± ^) 




B 


TO = 2 


n = 1 


(2,1) 


A R B + a Ra 


n (3)„(i) 


AB 


to = n + 1 


n = 1 


(2,1) 



Table 2: 



independent variables 


lapse choice 


constant of motion 


(B, AB) 


(2,1): Taub 


A'B^G^^ 


(AB) 


(2,-1): Brill 


A( 4 G^ ± ) 


(A, AB) 


(0,-1): Misner 


A 2 ( 4 G^±_) 



Table 3: 

The same problem may be viewed directly from the Einstein equations. 42 
Defining N 2 = A m B n , g D = diag(A, A, B), 4 R\ = 4 R 2 2 = 4 R A , and 4 i? 3 3 = 
4 Rb for the locally rotationally symmetric case, one has 

2N 2 ( 4 G- L ±) = -(hi A)' Qn By - ^lnA)' 2 + N 2 A~ 2 ( \n^ 2 B - 2n^n^A) , 
2A^ 2 ( 4 i? A ) = (In A)" - (In A)' {In N A' 1 B- 1 / 2 )- + N 2 A~ 2 (-n (3 > 2 B + 2n^ A) , 
2N 2 ( 4 R B ) = (In BY - (In By (In N A' 1 B- 1 / 2 )- +N 2 A' 2 (n^ 2 B) . 

(4.6) 

Different choices of (to, n) lead to decoupling of different products of A and B. 
Taub's 32 choice (2, 1) decouples (AB) -1 and i?, Brill's 43 choice (2, -1) decouples 
A -1 / 2 and B, while Misner's 44 choice (0, —1) decouples AB and A, the special 
exponents yielding the simplest form for the solutions. However, the decoupling 
of the evolution equations requires adding appropriate multiples of the super- 
Hamiltonian constraint to the Einstein evolution equations. Table 2 summarizes 
the decoupling possibilities for the evolution equations for the semisimple case. 

The first column of Table 2 lists those linear combinations of the three inde- 
pendent spacetime curvatures 4 Ra, 4 Rb, and 4 G ± ± which contain only one of 
the two independent spatial curvatures R 1 i+R 3 % (<~ n^n^) and i? 3 3 (<~ rj( 3 ) 2 ). 
The second column indicates which spatial curvature appears while the third 
column indicates the variable whose natural log appears in the second deriva- 
tive term. The fourth column gives the condition that the quadratic derivative 
term only involve this variable and the fifth column gives the condition that 
the quadratic term be absent. The last column gives the values of (m, n) for 
satisfying both conditions. 

Table 3 describes the choice of variables for which both evolution equations 
decouple in the semisimple case; the form of the constant of the motion may 
be derived from the evolution equations. In the nonsemisimple case of Bianchi 
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type II (n^ 1 ) — 0, ^ 0), the vanishing of the n^n^ curvature term leads to 
the additional lapse choice (1,-1) of Bradley and Sviestens 45 ' 46 for which both 
evolution equations decouple. When = = 0, one obtains the abelian 
case of Bianchi type I and the solutions are equivalent to the Kasner solution 
expressed with different time functions. 

Why give so much attention to a vacuum solution of d = 4 general relativity 
in an article devoted to higher dimensional cosmological models? The answer 
is that the Taub family of solutions are not content to remain in four spacetime 
dimensions but keep reappearing in higher dimensional models. 

For solutions of the Ricci evolution equations, the Taub time gauge Hamilto- 
nian H ~ g( 4 G ± ±) is a constant function of the parameters of those solutions. 
For negative values of H, one must add something positive to H to obtain zero. 
This corresponds to a source with E a b = (= E^ a ) but gT^± = H source > 0. 
Such a source is a homogeneous stiff perfect fluid moving orthogonally to the 
space sections, which in turn is equivalent to a homogeneous massless scalar 
field. The latter is equivalent by conformal transformations to a solution of the 
Brans-Dicke field equations. 9 In other words, by "varying the parameters" 42 
which occur in the solution of the Ricci evolution equations away from the vac- 
uum supcr-Hamiltonian constraint values, one obtains new solutions which may 
be interpreted as containing either or both of these equivalent sources or as a 
solution of the Brans-Dicke theory. (Such a variation of parameters applied to 
other choices of evolution equations leads to the introduction of a locally rota- 
tionally symmetric electromagnetic field source and the Brill generalization of 
the Taub solutions, 43 to which stiff perfect fluids or scalar fields may be added, 42 
or a cosmological constant. 47 ) 

These are still d = 4 solutions. However, the massless scalar field source 
or Brans-Dicke theory is equivalent to a d = 5 Kaluza-Klein theory with an 
additional flat dimension. 9 More scalar fields lead to higher dimensional Kaluza- 
Klein models, always with flat extra dimensions. 

In these examples, essentially only gravitational degrees of freedom have been 
considered. As an example of a model with nongravitational degrees of freedom, 
consider d = 11 supcrgravity with the Freund-Rubin ansatz. 10,11 Although al- 
ready out of fashion, it is a cute example of another appearance of the Taub 
solutions. Spatially homogeneous solutions of the bosonic sector of the theory 
have been considered by Freund, 10 Demaret et al 20 ' 21 and Lorentz-Petzold. 22 
The spacetime is assumed to be of the form M4 x Mr, where M4 <~ R x M3 
and the spatial metric is decomposable, i.e., equivalent to independent metrics 
on the factor manifolds M3 and M7 which are orthogonal in the product. The 
spatial metric matrix g is then in block diagonal form with matrix blocks g 3 
and g 7 , each of which may be decomposed into scale invariant and scale factor 
parts 

g 3 <-> (93, g 3 ) , §7 <"> (97, £7) • ( 4 -7) 

The Freund-Rubin ansatz which solves the field equations for the 3-form 
A of the theory assumes that the 4-form F — dA has as its only nonvanishing 



12 



components (let a, b = 1,2,3) 

F°^~ c = f D g, 1/2 , Wrf« l ~ c = -{Ngl' 2 )- 1 ' 2 ^ , e 0123 = 1 , 

(4.8) 

so that 

F^-F ^ = - \F a p 1& F<*W , (4.9) 

where fn is a constant (as in Dcmaret et al 20 ). 

This leaves an anisotropic source for the Einstein equations. The energy- 
momentum tensor or its Ricci equivalent form 

T a p = |;8(8F 2 Q/ 3 - g a pF 2 ) , E a(i = \{F 2 a[i - j2g a pF ) , 

F 2 aP = F ai& eF^ , F 2 ee F af376 F a ^ s , 

reduces to 

T% = -2/ D 2 37 - 1 , T\ = -2f D 2 g 7 - 1 5\ , T\ = 2f D 2 g 7 -H\ , 

(4.11) 

while the source super-Hamiltonian is 

jjsource = Af D 2 {g Z /g 7 ) 1 ' 2 , _ dlJ source = gl/^ab^ ^ ^ 

Note that by replacing F by its expression in terms of the metric leads to an 
effective potential which does not necessarily yield the correct equivalent energy- 
momentum upon variation. For the Freund-Rubin source no problem arises, but 
for a Freund-Rubin-Englcrt 12 source, one must introduce a nonpotential force 
to compensate. 

In Taub time gauge the source potential is 

N[J source = ( 2 . 10 ■ 9)g 1 / 2 U' mrCe = ( 2 ■ 10 ■ 9) ■ 4f D 2 g 3 . (4.13) 

Thus only 33 is affected by the source, the variables orthogonal to 173 continuing 
to satisfy the vacuum equations. To discuss the evolution of the metric one 
must make assumptions about the homogeneity groups of M3 and M 7 . If the 
spatial curvature on M3 is such that the equation for 53 is linear, i.e., the Taub 
time gauge potential is linear in g 3 , then adding the source only changes the 
coefficient, i.e., "varies the parameters" in the evolution equations, which have 
the same solutions but a different dependence on the parameters. 

In the simplest case M 3 and M 7 are both flat and the vacuum dynamics is 
just "free motion." With the Freund-Rubin source, the variables orthogonal to 53 
remain free and the single degree on freedom g% decouples from them leading to 
the 1-dimensional exponential scattering problem (in In #3); the solution is the 
Bianchi type II Taub solution involving one nontrivial curvature term, apart 
from a variation of parameters. The variables 53 and g 7 are not orthogonal, 
however. One must introduce 

g = e 2(/3 3 e 3 +/3 3 - L e 3i ) f«3 0^ e 3 = diag(l, 1, 1, -±, . . . , -±) 



, 3 = , g 7 = e'(-i^) 



g 7 y ' e 3± =diag(0, 0,0, !,...,!) ' (4U) 
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The (3 3 degree of freedom decouples from f3 3± which together with the remaining 
anisotropy variables remains free. The picture in the f3 3 -f3 3± plane is exactly 
the Bianchi type II Taub solution with a variation of parameters. In this plane, 
ln<?7 has slope ^; one sees that free motion in 31 - (parallel to the vertical 
potential contours) leads to monotonic expansion or contraction in 177. The 
initial asymptotic free state scatters off the static potential into a final free 
state with a velocity reflected from the potential contour direction, leading to 
an expansion from zero to a maximum and then recollapse in g 3 . Although 
this model 20 is not realistic, it does exhibit phases of preferential expansion or 
contraction of the ordinary and extra dimensions. (The solutions of the Ricci 
evolution equations for g 3 in Taub time gauge are in fact the same as the d = 4 
spatially flat case with a positive cosmological constant, but the proper time 
is a different function of the Taub time in the two cases due to the additional 
dimensions.) 

If instead one gives M 3 isotropic spatial curvature by letting its metric be the 
diagonal Bianchi type V metric, some special solutions may be obtained from 
the d — 4 vacuum solution (found by Joseph 48 ) by a variation of parameters. 
This is easily seen by considering this solution in Taub time gauge N 2 ~ g 3 
instead of the usual time gauge N 2 = (733 in which that solution is usually 
presented (discussed for spatially homogeneous models by Siklos 49 ). The free 
motion in the latter time gauge is represented by hyperbolic tangents instead 
of exponentials as occur in Taub time gauge, while the single curvature driven 
variable is described by hyperbolic sines in both gauges for this case. In Taub 
time gauge one sees that the Joseph solution is again the Taub solution with 
one curvature term (Bianchi type II) in disguised form. 

For a diagonal Bianchi type V metric (with C a jj- = a<5|") one has 

g 3 = e 2/3° e 2(/3+e ++ /3-e_) f e+ =diag(l,l,-2) , e_ = V3 diag(l, -1, 0) , 

(4.15) 

but the supermomentum constraint suppresses the (3 + degree of freedom (one 
may set /3 + = 0), resulting in a Hamiltonian system for the remaining two 
variables with g 33 = g 3 1 / 3 . In Taub time gauge the curvature potential is 

Nu curv = (2 . 1Q . 9 ) 57ff35 33 a 2 , = (2 . 1Q . 9) 5753 2/ V • (4.16) 

If one imposes that gj be proportional to 53 1 ^ 3 , this potential becomes pro- 
portional to #3 and simply adds on to the Freund-Rubin potential. The only 
difference from the vacuum case is that the coefficient of g 3 in this potential is 
different; the same general solution is obtained with a variation of parameters, 
i.e., the Taub solution occurs again. Apparently from the work of Dcmaret 
et al 20 and Lorentz-Petzold 22 this additional condition is compatible with the 
dynamics. If one could consistently impose the constraint that this potential de- 
pend on j3 3± only, then one would obtain a system equivalent to the semisimple 
Taub solution which describes the evolution in the presence of two exponential 
potentials depending on orthogonal variables. In general the two potentials are 
not orthogonal and the evolution equations do not decouple in Taub time gauge; 
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perhaps a direct analysis of the decoupling possibilities with power law lapses 
would lead to the solution. 

The Joseph solution is itself a special case of the Ellis-MacCallum diagonal 
type Vlh solution, 50 obtained by Lie algebra contraction of that solution. The 
same variation of parameters extends to this case and the whole family of d = 
4 Taub solutions, as well as the type VIo analog of the locally rotationally 
symmetric Taub metrics and the Kantowski-Sachs geometry, whose vacuum 
solution is related by analytic continuation to the Bianchi type III (= VI_i) 
value of the Ellis-MacCallum solution. These are discussed by Lorentz. 22 All of 
these are disguised versions of the Taub family of solutions which represents the 
dynamics of one or two nontrivial gravitational degrees of freedom, independent 
of the dimension. This variation of parameter idea is not special to d = 11 
supergravity in the higher dimensional context, but as already mentioned above, 
has implications for other theories as well. 

Consider the general case of a decomposable spatial metric from the point 
of view of the Ricci evolution equations 

«-(? L)- R2 \)' <«') 

Here Ri and R2 are the Ricci tensor mixed component matrices of the individual 
homogeneous metrics on the factor manifolds Mi and M 2 ; the source energy- 
momentum must also be in block diagonal form for consistency with the Einstein 
equations. The Ricci evolution equations are 

= a;- 1 (a;g- 1 g)-+2Ar 2 (R-E) , x = N~ 1 g 1 ^ 2 . (4.18) 

and in Taub time gauge x ~ const 

0=(g- 1 g )-+(2-10-9) 2 -2.g 1 .g 2 (R-E) . (4.19) 

In the diagonal case g = g D , the time derivative term is simply (lng^)". 

The only coupling between g 1 and g 2 in this time gauge occurs through the 
factor </i<72 in the driving term. When R 2 — E 2 = 0, for example, g 2 has free 
dynamics and one can rescale gi by a power of <? 2 to decouple the remaining 
evolution equations from g 2 

§1=52%, C=(£>i-l) _1 , I>i= dim Mi, 

o — — ( 4 -20) 
= (ir 1 gi)'+(2-10-9) 2 -25 1 (R 1 -E 1 ) . 

This redefinition of variables leads to the Taub time Ricci evolution equations 
for the spacetime R x Mi with spatial metric matrix g 1 . The conserved energy 
associated with those equations, however, must be nonzero due to the contribu- 
tion of the free motion in the extra dimensions to the higher dimensional energy. 
This was noted by Barrow and Stein-Schabes for a d — 5 Kaluza-Klein gener- 
alization of the Bianchi type IX models, for which generalized Taub solutions 
were mentioned. 23 
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When the spatial metric is not decomposable, but still has the block di- 
agonal form (4.17), a number of decoupling possibilities occur depending on 
the form of the driving term. Suppose that g 2 is the 1-dimensional matrix 
(9dd) so 52 = gDD- Then the gauge N = goD 1 ^ 2 leads to x = gi 1 ^ 2 ; if 
iV 2 (R— E) is independent of gDD, then g 1 immediately decouples from gDD, 
while In gDD satisfies a linear equation with a g x dependent source. When the 
trace of 7V 2 (Ri — Ei) depends only on x, taking the trace of the first block of the 
Ricci evolution equations leads to a decoupled equation for x, the kinetic term 
reducing to 2x~ 1 x. If N 2 Tr(Ri — Ei) is a constant, say — x 2 , this equation is 
in fact x — y^x = 0, leading to exponential or hyperbolic solutions when x 2 > 0, 
trigonometric solutions when % 2 < and linear solutions when x = 0. (This is 
again the 1-dimensional scattering problem, this time in the variable lnx.) This 
situation occurs, for example, if the spatial section is the semidirect product 
of a (D — l)-dimcnsional abclian subgroup and a 1-dimensional group of auto- 
morphisms of that subgroup, the 1-dimensional group being associated with the 
Dth direction. The block diagonal ("symmetric case") D = 3 nonsemisimple 
models are of this type, and this feature allows a uniform treatment in the time 
gauge N = 533 1 / 2 of most of the known vacuum solutions of this class, which 
are related by variation of parameters to scalar field (and hence higher dimen- 
sional Kaluza-Klcin), stiff perfect fluid (including tilt along the distinguished 
direction) and spatially self-similar generalizations. 42 Some solutions of the evo- 
lution equations, not allowed by the original supcr-Hamiltonian constraint, may 
satisfy that constraint after a variation of parameters, leading to solutions with 
no analog in the original case. Such solutions occur when a certain parameter 
changes sign. 

For x to decouple, it is sufficient that 7V 2 Tr(Ri — Ei) depend only on x. 
This occurs if this term arises from a curvature term associated with an isotropic 
coset space factor manifold to which the Z?th direction belongs. Fricdmann- 
Robcrtson- Walker or Kantowski-Sachs factors or even semisimple Taub factors 
lead to this behavior. In the D = 3 semisimple Taub case, the group manifold S3, 
as an Si fiber bundle over S2, 51 is locally a product, hence choosing N = A 1 / 2 
(i.e., (m, n) = (1,0)), where A is associated with the isotropic 5*2 base of the fiber 
bundle, leads to a decoupled equation for x = (AB) 1 / 2 . One may also rescale 
the lapse by a power of x without loosing this decoupling (i.e., m = n + 1). For 
example, in the Taub case, the lapse N = x~ 2 A = B~ x leads to the (0,-1) 
decoupling of x, yielding quadratic solutions. 

These ideas were used by Lorentz-Petzold (without being fully understood) 
in his discussion of some supergravity solutions 22 and Brans-Dicke solutions. 52 
The latter are equivalent to a d = 5 Kaluza-Klcin model in the vacuum case, 
and his choice of decoupling variable "5" involving the scalar field is exactly 
the variable x in the higher dimensional formulation. Adding perfect fluids 
with certain equations of state leads to a variation of parameters in the evo- 
lution equations without breaking this decoupling. All exact solutions can be 
understood in terms of the variation of parameters idea. 

Another situation arises when the lapse depends only on g, say N = f^ 1 ^ ^ 2 
so x = I 2 . The trace of the Ricci evolution equations then gives a decoupled 
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equation for x as long as N 2 Tr(R — E) is a constant, say — x 2 ; the equation 
is then x — ( 2 x 2x — 0, again having the same solutions as discussed above 
depending on the signs of ( 2 and x 2 - For example, in the d = 4 Bianchi 
type I perfect fluid with equation of state p — (7— l)p has E <~ <?~ 7 / 2 and hence 
N <~ g 7 / 4 makes iV 2 E constant; a cosmological constant alone has E constant, 
so in this latter case x has the same solutions as above but in proper time 
gauge. On the other hand for the Bianchi type I fluid case, one can also choose 
N = (^t^ 1 )/ 2 so that x ~ g 1 ' 1 ^ 2 and 7V 2 E ~ which leads to quadratic 
solutions for x; this was an intermediate step in the work of Jacobs. 53 

All these exact solutions are characterized by the fact that the source energy- 
momentum (including a possible cosmological constant) may be represented in 
terms of the metric and constants of the motion, leading to an entirely geomet- 
ric system. This occurs only when the source has either discrete or continuous 
additional symmetry. 37 In this class of models at most two "nontrivial" gravi- 
tational modes can be excited, as in the semisimple Taub case, the existence of 
exact solutions depending at least on a partial decoupling of the modes. 

5 Concluding Remarks on Exact Solutions 

At this point the overview of existing special exact solutions has shown us that 
they all depend crucially on decoupling ideas that involve the DeWitt geometry, 
the symmetry of the potentials or driving forces, the choice of time variable and 
the choice of evolution equations. The variation of parameters idea and its rela- 
tion to the constraints turns out to be a very important one, sometimes involving 
the parameters of the solutions of the equations and sometimes the parameters 
of the equations themselves. Although the super-Hamiltonian constraint has 
been emphasized because of the models discussed here, the supermomentum 
constraint plays a similar role when models are considered which excite the 
supermomentum. 42 All of these ideas seem never to be systematically applied. 
Most special solutions are found almost by accident and never later understood 
in a broader context. Studying the appropriate sections of the Exact Solutions 
book 54 with these ideas in mind, one begins to see some order in the seemingly 
unrelated list of exact solutions in ordinary spacetimc. In higher dimensions the 
situation is similar. 

But are particular exact cosmological solutions really significant? What do 
we gain by finding new very special solutions of some particular set of field 
equations? This is a timely question, given that the exact solutions industry 
often seems to be only weakly coupled to reality. One might imagine that 
special solutions might be interesting if they exhibit some "stable" property 
of the field equations in some sense or if the solutions exhibit some "typical" 
behavior during certain phases of the evolution of the universe, at least within 
the restricted symmetry class considered if not in a larger context. 

For the spatially homogeneous case, ordinary differential equations permit a 
sophisticated analysis of the qualitative behavior of their solutions. 55-58 In an 
appropriate formulation of the field equations, certain special exact solutions 
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play the role of "singular points" and "separatrixes" between singular points on 
an appropriate phase space. In the context of the parameter space of all symme- 
try types of a given dimension, exact solutions of a lesser symmetry type often 
characterize phases of evolution of a higher symmetry type. 38 Singular points are 
associated with "exact power law" solutions 59,60 of the field equations, which in 
turn are geometrically characterized by self-similar evolution, i.e., the existence 
of additional spacetimc symmetry in the form of a homothetic symmetry which 
shifts the homogeneous hypersurfaces along the time direction. This reduces the 
field equations to algebraic equations in certain parameters which determine the 
metric completely The simplest of these is the Kasner solution. The Taub so- 
lutions and their generalizations act as separatrixes in this picture. These kinds 
of solutions are important in understanding typical behavior at very early times 
toward the big bang and at very late times coming out of the big bang. 

Perhaps some of the energy devoted to finding particular exact solutions 
(a tree level activity) should be shifted to understanding some of the broader 
questions that arise in this area of research. The average level of "new exact 
solution" papers seems to have remained close to the level of two decades ago 
when this game first drew widespread participation. As we have seen, those 
exact solutions which exist are very special and have very nice mathematical 
properties which relate them to one another due to the richness of the Einstein 
equations. If these simple ideas are not commonly understood, how can one 
hope to have reasonable progress on more complicated questions? 

6 Chaos? 

One deeper mathematical question which seems to fascinate many people is the 
extent to which the "oscillatory approach to the initial singularity" found in 
the d = 4 spatially homogeneous vacuum semisimple case and in the (nonstiff ) 
perfect fluid case with anisotropic spatial curvature might be a "general prop- 
erty" of cosmological solutions of the classical Einstein equations. 28,29 This has 
aroused considerable controversy, more recently acquiring the trendier name of 
gravitational chaos 29 and spreading into higher dimensional purely gravitational 
models. 23 ^ 27,30,31 

In the d = 4 semisimple case, the vacuum supermomentum constraints re- 
strict the spatial metric to be diagonalizablc, the three spatial gauge degrees 
of freedom being intimately connected with the remaining three "offdiagonal 
modes" which are not excited in vacuum. (Diagonalization depends on the 
choice of the Lie algebra basis; the Killing metric must be diagonal in a basis 
in which the spatial metric is diagonal. Such a frame is then also an eigenbasis 
of the extrinsic curvature, namely a "Kasner frame" 28,46 .) In Taub time gauge 
the spatial curvature potentials at fixed a have "essentially closed" potential 
contours on the flat subspace of the two (3 variables, and for increasing values of 
the potential, these contour lines resemble more and more closely straight line 
contours joined together at vertices which are either closed or have open chan- 
nels which run out to infinity with a width which goes to zero. The "straight" 
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contours are exponential and move outward with a-speed ^ as a — > 0, while 
the system point moves with unit a-speed (a null geodesic) when the potential 
is negligible compared to the kinetic energy. Thus the system point continually 
overtakes and scatters off the receding exponential potentials as long as the rel- 
ative velocity between the system point in its free state and some exponential 
potential is always negative (approach). In fact the geometry is such that the 
relative velocity between the system point and the exponential potentials which 
is is chasing is nonpositivc in every direction but the vertex directions, where 
it vanishes. Motion exactly along a vertex direction only occurs for the Taub 
solutions which describe motion down the center of an open channel extending 
to infinity. Thus in general one has the Mixmaster behavior in which the sys- 
tem point rattles around indefinitely within the expanding potential as a — > 0. 
On the other hand the nonscmisimple vacuum models have an open set of di- 
rections in which the system point may escape, i.e., the absence of at least one 
of the "straight" potential terms leads to an open set of directions in which the 
potential asymptotically goes to zero, and the system point eventually ends up 
in an asymptotic free state. 

In the group manifold case, which is more general than the coset space case, 
each "straight" potential contour line is associated with the structure constant 
component C a bc with a ^ b ^ c, which enters as a quadratic factor in the 
potential. 1 For the d = 4 semisimple case, all such terms are necessarily present, 
leading to the essentially closed potential contours. In higher dimensions several 
crucial differences occur. One still needs semisimplicity for the existence of 
essentially closed potential contours; these potentials are best described in a 
frame in which the Killing metric is diagonal as in the d = 4 case. However, 
the Jacobi identity no longer allows all structure constants C a b c with distinct 
indices to be nonzero, removing certain potential terms. 

Furthermore, only some of the offdiagonal modes are restricted by the su- 
permomentum constraints, so the vacuum case in no longer diagonalizable and 
certain offdiagonal modes are nontrivial in the sense that the curvature potential 
necessarily depends on them and they cannot be suppressed without additional 
symmetry imposition, much like that which reduces the Mixmaster models to 
the essentially different Taub models. This leads to the excitation of nonzero 
curvature modes of the (rescaled) DeWitt metric. For those models which are 
diagonalizable and therefore no longer general, the absence of all possible struc- 
ture constants with distinct indices apparently leads to a vertex geometry such 
that an open set of vertex directions exists along which the "free" system point 
can chase the potential with positive relative velocity (recession); the system 
point then eventually ends up in an asymptotic free state. 26 ' 27 

Without imposing the Jacobi identity, Demaret et al 30 have shown that this 
happens only for d > 10, relevant to the general inhomogeneous case which still 
remains somewhat controversial. For the homogeneous case, the Jacobi iden- 
tities seem to reduce the critical dimension to d > 4 in diagonalizable models. 
However, in nondiagonal models, the eigenvectors of the extrinsic curvature no 
longer coincide with the time independent invariant frame vectors and the trans- 
formation to the eigenvector frame induces more nonzero structure functions. 
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Since it is the eigenvector frame (Kasner frame) structure functions which are 
relevant to the existence of an asymptotic free state, the question still remains 
open for those dimensions d € [5, 10] for which semisimple groups exist. 

Of course chaos even in d = 4 dimensions is sabotaged by the presence of a 
scalar field since it contributes a term to the super-Hamiltonian which reduces 
the free motion speed of the system point (timelike motion), therefore allowing 
an open set of vertex directions (or eventually all directions) which permit an 
asymptotic free state. 9 ' 37 Moreover, since the free motion phases are translations 
in a <~ Ing, one very quickly approaches the initial singularity g = 0, slamming 
into the Planck scale where the classical Einstein equations break down and the 
indefinite approach to the classical singularity is no longer particularly relevant. 
The classical chaos does lead to quantum consequences, however, at least at the 
naive level which allows calculation. 61 

7 Conclusion 

Whether one is interested in mathematical questions or more physical ques- 
tions, a broader point of view is extremely valuable. This article has focused 
on certain mathematical questions as a means of illustrating this point. Staring 
at equations or their solutions in a 2-dimensional representation is not partic- 
ularly illuminating. The rich geometrical structure of the Einstein equations 
allows one instead to exploit "hidden" relationships and visualize many con- 
cepts which remain somewhat obscure in their brute form representation. Even 
when visualization may be difficult or not particularly relevant, a perspective 
which views a given problem "from above" is certainly more effective in not 
only resolving that problem, but in understanding the result and placing it into 
context. The latter is, after all, more important than the former. 
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